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by judiciously choosing the feedback gains and the estimates of the nonlinear 
dynamics. In the first part of this paper, kinematics and dynamics modeling 
of the airships is presented. Euler angles and parameters are used 


in the formulation of this model and the technique of Computed Torque 

control is introduced. In the second part of the paper, we develop 
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1. INTRODUCTION 

Robust control of the autonomous airship has been the subject of considerable research over the last 
decade. The advantage of robust control over other control schemes is that robust control can be used 
in the presence of nonlinear dynamics. In the tracking control of the airship, the computed torque control law 
is commonly used. The computed torque control law has been shown to make the error system 
asymptotically stable, if the dynamics of the airship are exactly known. The Autonomous airship as typical 
lighter than air vehicle represents a unique and promising platform or various applications (civil and 
military). Such as telecommunications, broadcasting, the environmental protection (measurement of the air 
pollution, monitoring of the forests), scientific exploration [1]. There is a lot of key technologies to be 
realized for autonomous airship. One of the most important techniques is flight control system. Some 
contributions on dynamic modeling and control system have been reported in recent years [1]. The airships 
have a very low cost to maintain themselves in flight and have the possibility to fly in low altitude and with 
low speed. Airship is a member of family of under-actuated systems as it has fewer inputs than degrees of 
freedom and the principle of Archimedes applies in the air as well as under water. 

In this paper, the complete model (kinematic and dynamic) of the autonomous airship is presented, 
based on the Newton-Euler approach with analysis of gravity, buoyancy, aerodynamics, damping and thrust. 
This model is developing by [2]. The main objective of this article is the synthesis of stabilizing control laws 
in terms of translation and orientation for the airship. It’s a vehicle has the difficulty of the control because its 
complex, non-linear dynamics, multivariable, of coupled nature, particularly in its operation and sensitivity to 
the external disturbance. The backstepping controller is proposed to solve this problem (control). The control 
strategy is based on the decomposition of the original system into two subsystems: the first concerns 
the position control while the second is to control orientation and the velocity. 

There are many published works on the dynamic model and the control of lighter than air vehicles. 
In [3-5] motion is referenced to a system of orthogonal body axes located in the airship. This model was 
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formulated originally for a buoyant underwater vehicle [3]. It was modified later to take into account 
the specificity of the airship [4, 5]. This model was modified by [6] and the modification has the particularity 
that origin of the airship fixed frame is located in the centre of gravity. Recall that the centre of buoyancy is 
the centre of the airship volume. But in these papers, new dynamic model developed by [2] we used and we 
propose to control of this model. 

Many works are published in the control of the lighter than air vehicles, for example: 
in [7] the authors are proposed the linear control theory applied for the linearization. In [1] a fault-tolerant 
Backstepping control law is proposed. In [8] the authors proposed nonlinear adaptive controller design to 
stabilize the error system by Lyapunov direct method and Matrosov theorem. While in the cited works, they 
have used the decoupled linear models comprised of the longitudinal and lateral models. In [9] the authors 
have addressed the problem of designing tracking feedback control. In [10] they designed an adaptive 
feedback linearization control. In [11] an improved sliding mode control is presented. In [12] the dynamics 
in the vertical plane are analyzed and controlled, using maximal feedback linearization. In [13] a fuzzy 
adaptive backstepping control is presented, and the trajectory control method for under actuated stratospheric 
airship is proposed in [14] based on TLC theory. In [15] a backstepping controller was designed for attitude 
and velocity. In [16] a robust backstepping controller considering disturbance was proposed. In [17] a robust 
feedback controller is presented. In this paper, the computed torque controller and motion planning are 
combined to stabilize the Airship by using the point to point steering stabilization. Modeling is briefly 
described in Section 2. Section 3, describes computed torque controllers. In Section 4, simulation results are 
presented followed by robustness consideration in Section 5. Finally, some conclusions and future work are 
given in the last sectiona. 


2. MODELING OF THE AIRSHIP 

The airship Figure 1 studied here has a traditional ellipsoidal envelope. It aims to fly 
in the stratosphere. Buoyancy is provided by helium contained in the envelope. The avionics system, power 
system and payloads are equipped in a gondola fixed below the envelope. The aerodynamic control surfaces 
like rudders and elevators attach to the empennage surfaces. In this research, the up and down rudders move 
together, while the left and right elevators are differential and move separately. Therefore, the deflections of 
the rudders control the yaw movement whereas the elevators influence the pitch and rotation. The propellers 
are fixed on both sides of the gondola and provide main propulsive force for flight. 





Figure 1. Structure of the airship 


The dynamic modeling of lighter than air (airship) is deduced from Newton’s laws of classical 
mechanics. We are analyzing the motion of an airship with six degrees of freedom, whereas. We are 
considering the velocity of the wind is null. We described two references frames (inertial and body-fixed), 
then, the three-dimensional translational kinematic equations of an aerial vehicle taking into account the wind 
effect can be expressed as following (1): 


dx =V cos(w )cos(y) (1) 


dy =V cos(yv)sin(y) 
dy =V sin(y) 
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Where X and Y define the horizontal position of the vehicle and Z denotes its altitude. The translation 
equations are followed as (2): 


Yaa | 
—(b-—mg)cos(7) 
(m+m,,)V (2) 
(L+Tsin(a))sin(c) 
= (m+m,,)Vcos(y) 
d1 Tcos(a)- D 
K m+m,, Panne H 


These equations are a result of aerodynamics, gravity, buoyancy, lift, drag and thrust forces 
moments, and are controlled by the thrust T ,the angle of attack a and the bank angle o [2]. Aero dynamical 
Force: This force includes the drag force D, opposite to V, and the lift force L, orthogonal to V. These forces 
are expressed by the following (3): 


L= c, (M,a)A,pV° 


(3) 
D= TC, (M,a) A,pV? 


b=p*Vol*g 


The lift coefficient is a linear function of the angle of attack whilst the drag coefficient 1s a quadratic 


function of © (M9) (4) [18]: 


C, (M,a)=C,,(M)+K,, (M)a 
Cp (M,a)=C,, +KC; 


(4) 


These equations have an important place in aerospace vehicle study, because they can be formulated 
from trimmed aerodynamic data and simple autopilot designs. Nevertheless, they give a realistic picture of 
the translational and rotational dynamics unless large angles and cross-coupling effects dominate 
the simulations. Trajectory studies, navigation, and guidance evaluations can be successfully executed with 
simulations of these equations. The limitations on thrust and velocity will be used in the path 
planning [18-20]. 


3. NONLINEAR COMPUTED TORQUE CONTROL 
3.1. For the dynamic model 

The velocity, the pitch and the yaw can be controlled by a feed forward controller Figure 2. 
From the model (2) one can see that through the velocity V depends on T and the other hand o and a are 
controlled through the angles Y, y respectively. 





ackstepping | 
controller | 






Pr 





: Backstepping a Y 


Figure 2. Computed torque controller for dynamic model 
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The computed torque algorithm to develop the control allowing the system to follow the desired 
trajectory, we consider the tracking error (6): 


e,=V-V (6) 


é,=V-V, (7) 


The nonlinear feedback control law that guarantees tracking of desired trajectory. Selecting Proportional plus 
derivative (PD) feedback for U (t) results in the PD-computed torque controller [3], the final control is (8): 


Tasin(a)+D 
—7(b-—mg)coscos(y) 
E +(V, -Ke -K„ė,)(m+m,) 


r vl~v (8) 
COS COS Q 
and the resulting linear error dynamics are (9): 
(Ñ, ~ Ke, E Kè, ) = 0 (9) 


According to linear system theory, convergence of the tracking error to zero is guaranteed. 
The yaw and pitch attitudes can be stabilized to a desired value with the following tracking feedback 
control (10-11). 


((L+ Tsin a)coso +(B-mg)cos 1)V 
+Vo(L +Tsin a.)sin o 

-(L+Tsin a) V coso 

-V}(B — mg)sin y 


+(¥, -K „e, -K,,é, }(m+m,;,)V" (10) 
T V cosa coso 


Q= 


V 
ies j (L+Tsina)sino 
—Vysiny 
(L+Tsina+Técos a) 
V sin o cos y 


V? cos? ¥ 


(L+Tsina)V cosocosy 


Aon rae —ky,ey — Kyzêy ) (11) 


= 


where Ky1, Ky2, Kyi, Ky2» Ky, and Ky; are the coefficients ensuring stability and are positive constants. 


3.2. For the kinematic model 
In Figure 3 is the motion through the axes X, Y and Z can be controlled by a feed forward controller. 
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Figure 3. Computed torque controller for kinematic model 


From the model (1) one can see that the motion through the axes Z depends on y, on the other hand 
V and Wcontrol the motion through the axes X,Y respectively. This leads to a Backstepping controller for 
which X, Y, Z are given by (12-14): 


X, —kyiex —kyyex 
+V (Ysin ¥ cosy +/sin y cos Y) (12) 
V= 
cos ¥ cos y 


le — ky ey E kyy | 


-(V cosy- Výsin y)sin ¥ 





ý- (13) 
V cos ¥ cos y 


. BA -kze —kz,€, —Vsin y] (14) 
y Z MMMM 
V cos y 


where Ky, Kx2, Ky1, Ky2, Kzįand Kz, are the coefficients ensuring stability and positive constants. 


4. SIMULATION 

To evaluate the designed control system, repetitive simulation tests were performed via numerical 
simulation. The control system was simulated using the variable step Runge—Kutta integrator in MATLAB. 
The airship is tested in simulation in order to validate motion planning algorithm considering the proposed 
Computed Torque controller. The general specifications of the vehicle are as follow: 

m = 18.375 (Kg); Vol = 15 (m? ); p = 1.225; 


M1, = 1.247 * p(Kg); m33 = 16.671 * p(Kg); 
m 

m22 = 17.219 * p(Kg); g = 9.8 (=); 

CDo = 0.5; CLo = 0.024; CL, = 0.937 

K = 1.4;A, = 2.84 


The parameters of the computed torque controller are selected as: 
For dynamic model: 


Kyi = 250, Ky2 = 130, Kw = 25, Kw = 2; 
Kyi = 16, Ky2 — 1.5 


For kinematic model: 
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Ky, = 16, Ky = 2.3, Ky, = 10.5, Ky, = 5.6, 
K714 = 15, Kz; = 6.9 


Figure 4 shows the simulation results of PD-CTC control with accurate parameters and without 
external disturbances and to confirm the robustness of proposed trajectory tracking control method. 
As illustrated in the above figures, the convergence to the trajectory is guaranteed after a transient behavior. 
It can be observed that the airship platform tends to the target point precisely, which demonstrates that 
the proposed approach succeeds in station keeping control for the airship platform. Figure 5 shows that 
the position errors and Figure 6 the translation dynamic (velocity and orientations) of the airship platform 
with accurate parameters and without external disturbances. The position errors are converging to zero, 
the translation dynamics converge to desired value, and we can see the change of the translations dynamics 
(angles and velocity) in the moment for this change in the trajectories (20s, 40s, 60s). The above mentioned 
results demonstrate that the PD-CTC control is accomplished with precision using the designed control 
system. In conclusion, the simulation results verify that the trajectory tracking control method designed 
in this paper for airship is efficient. 
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Figure 5. The position errors of the airship without disturbances 
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Figure 6. The translation dynamic (velocity and orientations) of the airship without disturbances 


CONCLUSION 
In this paper, we presented synthesis of stabilizing control laws by computed torque approach. 


Firstly, we start by the development of the dynamic model of the airship taking into account the different 
physics phenomena, using Newton Euler approach and for complex non-linear movement of blimps we can 
design its computed torque controllers. As prospects we hope to develop other control techniques in order to 
improve the performances and to implement them on a real system. 


REFERENCES 


[1] 
[2] 
[3] 
[4] 


[5] 
[6] 
[7] 
[8] 


Guoqiang Wu, BaoJun Lin, Shancong Zhang, “Fault-Tolerant Backstepping Attitude Control for Autonomous 
Airship with Sensor Failure,” Procedia Engineering, vol 29, p 2022 — 2027, 2012. 

Y. Bestaoui Sebbane, “Lighter than air robots intelligent systems control and automation,” journal of Guidance, 
Control, and Dynamics, vol. 32, no. 4, Jul-August 2009. 

L. Beji & A. Abichou, “Tracking control of trim trajectories of a blimp for ascent and descent flight maneuvers,” 
International Journal of Control, International Journal of Control, vol. 78, no 10, p. 706-719, 2005. 

Wu Yongmei, Zhu Ming, Zuo Zongyu, Zheng Zewei, “Trajectory Tracking of a High Altitude Unmanned Airship 
Based on Adaptive Feedback Linearization,” International Conference on Mechatronic Science, Electric 
Engineering and Computer, August 19-22, Jilin, China, 2011. 

Yueneng, Yang, Wei, Zheng, et Jie, Wu. "Sliding mode control for a near space autonomous airship, ” In: Electric 
Information and Control Engineering (ICEICE), International Conference on, IEEE, p 3576-3579, 2011. 

Xiaotao Wu, Claude H. Moog, Luis Alejandro Marquez-Martinez, Yueming Hu, “Full model of a buoyancy-driven 
airship and its control in the vertical plane,” Aerospace Science and Technology Elsevier, vol 26, p 138-152, 2013. 
Yueneng Yang, Jie Wu, Wei Zheng, “Station-keeping control for a stratospheric airship platform via fuzzy adaptive 
backstepping approach,” Advances in Space Research, vol51, p1157—1167, 2013. 

Zewei Zheng, Wei Huo, Zhe Wu, “Trajectory tracking control for Underactuated Stratospheric Airship,” Advances 
in Space Research, vol 50, p 906-917, 2012. 

Y. Zhang, W. Qu. “Stabilization and trajectory tracking of autonomous airship’s planar motion,” Journal of Systems 
Engineering and Electronics, vol 19, no 5, p 974-981, 2008. 

K. S. Kim and Y. Kim, “Robust backstepping control for slew maneuver using nonlinear tracking function,” IEEE 
Trans. Control Syst. Techno, vol. 11, no. 6, pp. 822—829, Nov. 2003. 

X. Wang, X. Shan, “Study of robust control for stratosphere airship attitude,” Journal of System Simulation, 
vol 5, p18-5, 2006. 

Bestaoui, Y. “Mission plan under uncertainty for an autonomous aircraft,’ Proceedings of the Institution of 
Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 224, no 12, p. 1297-1307, 2010. 

B. Siciliano and O. Khatib, Springer handbook of robotics: Springer-Verlag New York Inc, 2008. 

Farzin Piltan and all, “Design Artificial Nonlinear Robust Controller Based on CTLC and FSMC with Tunable 
Gain,” International Journal of Robotics and Automation (IJRA), Volume (2): Issue (3): 2011. 


Int J Rob & Autom, Vol. 8, No. 1, March 2019: 44-51 


Int J Rob & Autom ISSN: 2089-4856 g 51 


[15] 
[16] 


[17] 
[18] 
[19] 


[20] 


T. I. Fossen, Guidance and Control of Ocean Vehicles John Wiley Sons Ltd, Chichester, (1994). 

E. Hygounec, P. Soueres, S. Lacroix, “Modélisation d’un ballon dirigeable, Etude de la cinématique 
et dynamique,” CNRS report, 426, LAAS, Toulouse, France, October 2000. 

G. A. Khoury, J. D. Gillet, eds. “Airship technology,” Cambridge University press, (1999). 

L. Beji, A. Abichou, Y. Bestaoui, “Stabilization of a Nonlinear Underactuated Autonomous Airship-A Combined 
Averaging and Backstepping Approach,” third international Workshop on Robot Motion and Control, 
November 9-11, 2002. 

Y. Bestaoui, T. Hamel, “Dynamic Modeling of Small Autonomous Blimps,” Conference on Methods and models 
in automation and robotics, Miedzyzdroje, Poland, vol. 2, pp.579-584, 2000. 

Zhang Yan, Qu Wei-Dong, Xi Yu-Geng, Cai Zi-Li, “Adaptive Stabilization and Trajectory Tracking of Airship with 
Neutral Buoyancy,” Vol. 34, No. 11 November 2008. 


APPENDIX 

[x,y,z]! Position vector 

[W, y]" Euler angles yaw, pitch 

V Velocity of airship 

m Masse of the vehicle 

T Force developed by propulsion 

[o , a]! Angles of the rudders and elevators 
Ar Setting of the reference area 

p Density of the air 

M Mach number 

Vol Volume of the vehicle 

b = Force of buoyancy 

CL (M, a), Cp(M, a) Parameters of lift and drag 
L, D =Lift and drag forces 
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